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We devise a technique for defining and computing n-point functions in the context of a background- 
independent gravitational quantum field theory. We construct a tentative implementation of this 
technique in a perturbatively- finite loop/spinfoam model. 

The lack of a general technique for computing particle scattering amplitudes is a seriously missing ingredient in 
nonperturbative quantum gravity [1,2]. Various problems can be traced to this absence: the difficulty of deriving the 
low energy limit of a theory; of comparing alternative theories, such as alternative versions of the hamiltonian operator 
in loop quantum gravity (LQG) or different spinfoam models; or comparing the predictions of a theory with those 
of perturbative approaches to quantum gravity, such as perturbative string theory. Here we explore one possibility 
<^ • for defining a general formalism aimed at computing scattering amplitudes. We outline a calculation strategy, which 
can certainly be improved. Our interest is not in a particular theory, but rather in a general technology to be used 
C~\ . for analyzing different models. For concreteness, we implement this strategy in the context of a specific model, and 
present a well-defined and perturbatively finite expression, which, under substantial assumptions and approximations, 
might be used as a general covariant n-point function. 
^\ ' In conventional QFT, we can derive all scattering amplitudes from the n-point functions 

W{x u . . . , x n ) = Z- 1 ( D<p <f>{x x ) . . . cj>{x n ) e- iS M (1) 

\q ; J 

where the Xi, i — 1,. .. , n are points of the background spacetime, (j) is t ne quantum field, S[(j)] its action and Z 
is the integral of the sole exponential of the action. Alternatively, the n-point functions can be derived from their 
Euclidean continuations, defined by dropping the i factor in the above expression. The integral (1) is well-defined 
in perturbation theory or as a limit of a lattice regularization, under appropriate renormalization. A well-known 
difficulty of background independent quantum field theory is given by the fact that if we assume (1) to be well-defined 
| with general-covariant measure and action, then the n-point function is easily shown to be constant in spacetime (see 
H> for instance [3]). This is the difficulty we address here. 

5^ | Consider a spacetime region R such that the points Xi lie on its 3d boundary E. Call ip the restriction of the field 
(j> to E. Then (1) can be written in the form 

W(x 1 ,...,x n ) = Z- x J Di P i P {x 1 )...i P {x n ) W[<p^ W[<p,Y,} (2) 

where (in the euclidean case) 

W[tp,H] = f D<f> R e ~ s «[^] (3) 

is the integral over the fields in R bounded by the 3d field ip on E, and Sr is the restriction of the action to R; while 
W[ip, E] is the analogous quantity defined on the complement R of R in spacetime. The field boundary propagator 
(3) has been considered in [4] and [1] and studied in [5]. Assume that we are dealing with an interacting theory, 
approximated by a free (gaussian) theory S^ )^] in some regime, and that, within a certain approximation of the 
amplitude (1), the interaction term in the action can be restricted to R. Then we can replace W[^j, R] with its free 
theory equivalent 



Wbb.E] = / D0 R e- b RM ee ttsM. (4) 

This integral is gaussian and can be performed, giving a gaussian "boundary state" 'Fx; [ip] , determined by appropriate 
boundary conditions for the field at infinity. For instance, if we take R to be defined by t > 0, then $ f -o[ip] is the 
vacuum state in the functional Schrddinger representation. In general, we expect the boundary state to be given by 
some gaussian functional of the boundary field ip on E. 
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Consider a diffcomorphism invariant theory including the gravitational field. Assume that the equations above 
hold, in some appropriate sense. The field <j> represents the gravitational field, as well as any eventual matter field, 
and we assume action and measure to be diffeomorphism invariant. Two important facts follow [6]. First, because of 
diffcomorphism invariance the boundary propagator W[<p, S] is independent from (local deformations of) the surface 
S. Thus in gravity the left hand side of (3) reads W[^>]. Second, the geometry of the boundary surface S is not 
determined by a background geometry (there isn't any), but rather by the boundary gravitational field ip itself. 

We can obtain an indication on the possible forms of the boundary state in gravity from the free quantum theory 
of non-interacting gravitons on Minkowski space. If we take R to be t > 0, for instance, then ^t = o[y] must be 
approximated by the well-known Schrodinger vacuum wave functional of linearized gravity. This is a gaussian state 
picked around a classical geometry: the flat geometry of the t = surface in Minkowski space. In the case of a 
compact R, it is then reasonable to consider a gaussian boundary state ^f q [f] picked around some 3-geometry q of 
the boundary surface S. Thus, we may expect an expression of the form 



to approximate (1) when the interaction term can be neglected outside R. For this equation to be significant, we 
have to fix the meaning of the coordinates Xi, since the rest of the expression is generally covariant. There is an 
obvious choice: the points Xi can be defined with respect to the geometry q. For instance, if n = 4, = t® = and 
*3 = *2 — ^\ ( we use x = (A^O) we can t& ke q to be the geometry of a rectangular box of height T and side L and 
interpret Xi as proper distances from the boundaries of the box. In other words, the localization of the arguments of 
the n-point function can be defined with respect to geometry over which the boundary state is picked. Notice that Xi 
in (5) are then metric coordinates: they refer to gravitational field values. They are not anymore general-covariant 
coordinates as in (1). In this manner, we can give meaning to n-point functions in a background independent context. 

Physically, we can interpret J? as a finite spacetime region where a scattering experiment is performed. The 
quantities Xi are then relative distances and relative proper time separations, measured along the boundary of this 
region, and determined by (the mean value of) the gravitational field (hence the geometry) on this boundary. This is 
precisely the correct general- relativistic description of the position measurements in a realistic scattering experiment. 

In order to give (5) a fully well-defined meaning, and compute n-point functions concretely, we need four ingredients: 
(i) A proper definition of the space of the 3d fields ip integrated over, and a well-posed definition of the integration 
measure, (ii) An explicit expression for the boundary propagator VF[<p]. (iii) An explicit expression for the boundary 
state [<£>]. (iv) A definition of the field operator <p(x). In the following we analyze the status of these four ingredients 
in the loop and spinfoam approach to quantum gravity. We consider for simplicity pure gravity without matter. 

(i) In quantum theory, the boundary values of Feynman integrals can be taken to be the classical dynamical variables 
only if the corresponding operators have continuum spectrum. If the spectrum is discrete, the boundary values arc 
the quantum numbers that label a basis of eigenstates (see [1]). In our case, the boundary field ip represents the 
metric of a 3d surface. Let us assume here the results of LQG that the 3d metric is quantized [1]. Therefore we must 
replace the continuum gravitational field variable ip with the quantum numbers labelling a basis that diagonalizes 
some metric degrees of freedom. These can be taken to be (abstract) spin networks s, or s-knots. An s-knot is here 
an equivalence class under (extended [7]) diffeomorphisms of embedded spin networks S. An embedded spin network 
is a graph immersed in space, labeled with spins and intertwiners. The s-knots arc discrete [7]. Thus, we rewrite (5) 
in the form 



where the meaning of ip(x) will be specified later on. The discrete measure c(s) on the space of the s-knots is defined 
by the projection of the scalar product of the space of the embedded spin networks (c(s)=l) except for discrete 
symmetries of s.) For simplicity, and in order to match with the spinfoam formalism that we use below, we restrict 
here the space of the spin networks to the four-valent ones and we identify spin networks with the same graph, spins 
and intertwiners (i.e., we ignore knotting and linking). 

(ii) The boundary propagator W[s] is a now a function of a boundary spin network. A natural possibility is to 
identify it with the boundary propagator W[s] defined by the spin foam models [8]. For concreteness, let us choose 
here the model defined by the SO(A)/ SO(3) group field theory [9], which gives a perturbation expansion finite at all 
orders [10]. This is the model denoted GFT/C in [1]. The amplitude of a spin network s is given in this model by 




(5) 




(6) 



S 




(7) 
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Here $ is a function on [S'0(4)] 4 and the precise meaning of the (symbolic) integrals in the exponent is detailed in 
[8] and [1]. The quantity / s [$] is a polynomial in the field <&, determined by s. It is defined by picking one factor 



ai...ft4 



= jd 9l . ..dg 4 $( 3l , . . . ,.g 4 ) . . .R { it Wi {9A) ^...p A (8) 



per node of s, where v l is the intertwiner of the node and ji, ■ ■ ■ ,ji the colors of the adjacent links, and contracting 
the indices oti according to the connectivity of the graph of s. The expression (7) is well-defined and finite order 
by order in A. (The rigorous proof of this statement is complete up to certain degenerate graphs [8].) The explicit 
computation of W[s] is entirely combinatorial and can be performed in terms of combinations of n J Wigncr symbols 
[1]. For completeness, recall that the reason for the definition (7) is that the expansion of W[s) in A can be written 
as a sum over spinfoams bounded by the spin network s 

W[a] = J2 ( 9 ) 

8(7 = 8 

where the spinfoam amplitude A(a) is the Barrett-Crane discretization of the exponential of the Einstein- Hilbert action 
of the discrete four-geometry defined by the spinfoam a. Therefore the definition (7) of W[s] can be interpreted as a 
(background independent) discretization of the functional integral (3). 

(iii) An expression for the boundary state \& 9 [s] can be obtained from the analysis of the coherent states in LQG 
[11-13]. For concreteness, let us pick here Conrady's definition of a coherent state [13]. Other more refined expression 
could be used instead. Conrady has defined a state ^o[S] that describes the Minkoski vacuum as a function of 
embedded spin networks 5*, under certain approximations and assumptions. This function has the property of being 
picked on spin networks that are "weaves", namely that approximate a flat metric q when averaged over regions 
large compared to the Planck scale [14]. This vacuum state can be written as follows. Pick cartesian coordinates 
x a , a = 1, 2, 3, on a 3d surface equipped with a flat metric q and with total volume V. Fix a triangulation T of lattice 
spacing a, small compared to the Planck length l p in the metric q. Restrict the attention to embedded spin networks 
S living on T. Define the form factor of a spin network as 



F s b ^ = &T, E f dt f dt'e a {t)e\t')8{x-x v ), 



(10) 



where v are the vertices of the spin network S; x v their position; e : t n e°(t) the edges; e e v indicates the edges e 
adjacent to the vertex v; and e a — de a /dt. Its Fourier transform is F§ b (k) = V^ 1 J d 3 x c ~ lk ' x Fg b (x). Then 

* [S] =AAcxp[- -L £ \k\ \Ff(k)j e (j e + 1) - VV5 ab 6^\\ (11) 

where the momenta summed over are the discrete modes on the triangulation; j e is the spin associated to the edge 
e; TV is a normalization factor. To understand this construction, notice that if we consider the gravitational field 
associated to the spin network, (in the sense of the weaves) qs b (x) = Fg b (x) j e (j e + 1), then ^fo[S] = ^oiqs] where 

*o[<?] = / d3x f d3 y^ ab (S)-s ab )WA(S-y)(q ab (y)-s ab )] ^ 

is the Schodinger functional representation of the linearized vacuum state. Wa(x — y) is a lattice rcgularization of the 
vacuum covariance. We can extend this construction to a 3d (euclidean) rectangular boundary S simply taking the 
product of the Conrady states associated to each of the eight faces forming S. 

We need to carry this result over to the diffeomorphism invariant s-knot states. Given an abstract spin network s 
there will be in general one embedded spin network S(s) that maximizes the state *o['5]- We can then tentatively 
define ^>o[s] = ^q[S(s)]. Notice that if *f?o[S] is picked on weaves, then the diffeomorphism invariant state ^o[s] defined 
is picked on the corresponding ("weavy") discrete 3-geometries. The maximization condition can be interpreted as a 
gauge choice, picking the coordinate system in which the 3-geometry is closest to the euclidean metric. The gauge 
invariant state is then chosen to be the restriction of the state to this gauge surface. In the spirit of [13], we restrict 
to embedded spin networks S living on T. Given an s-knot s, there is only a discrete number of such spin networks 
that are in the class s: we choose S(s) that maximizes (11) among these. We expect this definition (possibly with 
an appropriate correction of the Conrady vacuum state) to converge for fine triangulations, making the background 
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structure chosen effectively irrelevant for a triangulation sufficiently finer than the Planck length. This construction 
provides a finite definition of 'I'oH, diffeomorphism invariant by definition. 

(iv) Finally, we need to define the field <f)(x) appearing in (6). Following [13] we write h^ b (x) = (qs\ s )( x ) — <5 ah )> 
where the point x is defined in terms of the boundary metric q and <Zg( s ) is defined above (12). An alternative, which 

we do not pursue here, is to derive h ab (x)(S) from the action of two SU(2) generators [15]. 

We can now bring together the various pieces discussed. To start with, consider a parallelepiped in 4d euclidcan 
space with hight T and cubic base of side L. Let £ be its boundary, equipped with the induced metric q. Fix a 
triangulation of X. The simplest choice is to start from a cubic triangulation of £, and to obtain a four-valent lattice, 
by splitting each (six-valent) vertex of the cubic lattice into two vertices. Replacing the various items discussed into 
the formal expression (5) we obtain 

W a ^- a ^(x u ...x n ;L,T) = ZzU2 c & K lbl (xi)---K nb "M W[s] (13) 

s 

where the sum is over all the s-knots that can be embedded in the triangulation. The normalization factor is the 
"vacuum to vacuum" amplitude Zt,t = ^2 S c(s) ty q [s]W[s\. (13) can be expanded in powers A", n is the number 
of vertices of the spinfoam, which is the number of 4-simplices in a simplicial complex dual to the spinfoam, if this 
exists. As a rough estimate, we can imagine each 4-simplex to have Planck size: if classical configurations dominate, 
the main contribution should come from n of the order of the 4-volume of the interaction region in Planck units. 

Hypotheses and approximations used to get to (13) are severe. But all quantities in (13) are well defined. The 
expression is probably finite at any order in A. We can thus take (13) as a tentative definition of an n-point function 
within the formalism of non-perturbative quantum gravity. More precisely, we can consider (13) as a tentative concrete 
definition of the quantity formally given by 

W aibl - a " b "(x u ...,x n ) = Z^ J Dg g aibl ( Xl )... <? a " 6 "(x„) e ~ SBH ^ (14) 

where Seh is the Einstein-Hilbert action, computed at relative spacetime distances x\, . . . , x n evaluated in terms of 
the mean value of the quantum gravitational field itself, on a box encircling the interaction region. 

The construction can probably be ameliorated and varied in a number of ways, and many issues remain open. 
There are important missing steps to get to the definition of quantities that can be interpreted as particle transition 
amplitudes. (On the physical interpretation of "particle" states defined on finite spacial regions, see [16].) The key 
questions are whether the expression (13) is indeed finite, convergent, and independent from the auxiliary structures 
uses to define it, when the triangulation is sufficiently finer than the Planck scale, and whether this construction leads, 
in a first approximation, to the general relativity scattering tree amplitudes. 
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